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Abstract 

Second order super symmetry transformations which involve a pair of complex con- 
jugate factorization energies and lead to real non-singular potentials are analyzed. 
The generation of complex potentials with real spectra is also studied. The the- 
ory is applied to the free particle, one-soliton well and one-dimensional harmonic 
oscillator. 
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1 Introduction 



The n-th order supersymmetric quantum mechanics (n-SUSY QM), which in- 
volves differential intertwining operators of order n, is a useful tool for generat- 
ing new solvable potentials [1,2,3,4]. Due to its simplicity, the 1-SUSY QM is 
the most explored; its nonsingular transformations produce partner potentials 
whose spectra can differ at most in the ground state energy level [5,6,7]. The 
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difficulty of 'modifying' the excited part of the spectrum has been surpassed 
through the 2-SUSY QM [8,9,10,11,12,13,14], which allows to 'create' two new 
levels €i, e 2 between two neighboring energies E iy E i+1 of the initial Hamilto- 
nian [15]. A similar treatment, implemented for periodic potentials [16], can 
be used to embed two bound states in a spectral gap above the lowest energy 
band [17,18]. In both situations the corresponding 2-SUSY transformations 
are irreducible, i.e., when obtained as the iteration of two 1-SUSY procedures 
they will involve always singular intermediate potentials. 

Here we will study a different set of 'irreducible' 2-SUSY transformations ap- 
plied to non-periodic potentials, which employs two complex conjugate factor- 
ization energies e±, e 2 , e 2 = e±. Irreducibility means now that the intermediate 
potential is complex although the final one is real. This problem has been ad- 
dressed previously [19], but up to our knowledge the conditions granting that 
the final potential will be regular have not been yet examined. We will show 
as well that the non-singular case leads to intermediate complex potentials 
having real spectra. These points constitute the subject of this letter, which 
has been organized as follows. In section 2 the second order SUSY transforma- 
tions with ei, e 2 € C, e 2 = ei will be analyzed. A prescription for avoiding the 
singularities in the new potential will be provided in section 3, while section 
4 will be devoted to study the intermediate complex potentials. Section 5 will 
deal with some particular examples as the free particle, one-soliton well and 
harmonic oscillator potential. 



2 Second order supersymmetric quantum mechanics 



The standard supersymmetric (SUSY) algebra with generators Qi, Q 2 (su- 
percharges) and H ss (SUSY Hamiltonian) reads: 



{Qj, Qk} = 5 jk H ss , [H ss , Qj] = 0, j, k = 1, 2, 



(1) 



where [•, •] denotes a commutator and {•, •} an anticommutator. The realization 

Qi = (<2 f + Q)/V2, Q 2 = {Q ] - Q)/(iV2) with 
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is called second order sypersymmetric quantum mechanics (2-SUSY QM). In 
this formalism H, H are two intertwined Schrodinger Hamiltonians: 



HA = AH, (5) 
and thus the real functions rj(x), j(x) are related with V(x), V(x) through: 



J-V = 2 V \ (7) 
(V- V)r] = 2V' + 2^' + V ", (8) 
(V - V)i = rjV + V" + 7" . (9) 

To decouple this system, substitute (7) in (8) and integrate to obtain 

>y = d-V + rf i /2-rf/2, (10) 

where d G K is a constant. By plugging (7,10) into (9), multiplying the result 
by 77 and performing the next integration one arrives to: 

W - (V) 2 / 2 + 2?7 2 ^ 2 /4 -ri'-V + dj+2c = 0, (11) 

c G K being another constant. This formalism is useful for generating new 
solvable potentials. To illustrate it, suppose that c, d are fixed and V(x) is 
an initial exactly solvable potential. The complete determination of V(x) in 
(7) requires to find solutions i](x) of the nonlinear second order differential 
equation (11). Let us find them using the Ansdtz [12]: 

rj'( x ) = r,\x) + 2(3{x)r ] {x) - 2^(x) , (12) 

where j3(x) and £(x) are to be determined. After substituting (12) into (11) we 
obtain a system of equations from which it follows that £ 2 = c. The essential 
part of the system is the Riccati equation: 

(3' + {3 2 = V-e, e = d + £. (13) 

As there exist two possible values of e, e\ — d + \fc and €2 = d — y/c (which 
coincide with the factorization energies in (3)), we indeed are dealing with two 
equations (13) whose solutions will be denoted by f3i(x), ^{x). This leads to 
a natural classification of the 2-SUSY transformations based on the sign of c, 
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which will be discussed elsewhere. Here, we restrict ourselves to the complex 
case for which c < and then e±, e 2 € C, e 2 = e±. Since V is real we can take 
Pi{x) = Pi(x), i.e., the problem reduces to solve the Riccati equation for 
In addition, using (12) we get two equivalent expressions for the real rj(x): 



?7 / = ?7 2 + 2/3 1 77-2iIm(ei), (14) 
ri' = ri 2 + 2(3 1 r] + 2ilm(e 1 ). (15) 

By subtracting both equations and solving for 77, we obtain: 



77 = Im(e 1 )/Im(/3 1 ) . (16) 
Once we know 77, the 2-SUSY partner potential V(x) is calculated using 

V = ^ + 2[Im(e 1 )/Im(/3 1 )]' . (17) 

Let us remark that the case we are dealing with has been previously explored 
[8]. However, we have not found any previous analysis on how to avoid the 
singularities in V(x), a phenomenon which seems almost unavoidable in the 
complex case [19]. 

3 The non-singular 2-SUSY potentials 



The simplest algorithms departing from and arriving at the exactly solvable 
potentials should avoid the singularities which might appear in V(x). Notice 
that a singular V(x) could be treated as a non-singular partner potential of 
V(x) in a restricted rc-domain. However, this would require at the end to solve 
the initial Schrodinger equation with modified boundary conditions loosing, 
in general, the solvability of H [20]. 

Let us rewrite first the formulae of section 2 in terms of solutions ui(x) of the 
Schrodinger equation arising from (13) by the change f3(x) = u'(x)/u(x) [14]: 



-u"(x) + V(x)u(x) = eu(x). (18) 

Hence r)(x) = — 2ilm(ei)\ui\ 2 /W(ui, Hi), where W(ui,U2) = u\v! 2 ~~ u 2u[ de- 
notes the Wronskian of u±, u 2 . From now on it is convenient to work with the 
real normalized Wronskian w(x) = W(ui,Ui)/[2ilm(ei)]. Therefore: 
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w \x) = \ Ul (x)\ 2 , (19) 
T)(x) = —w'(x)/w(x), (20) 
V( x ) = V(x) - 2[w'(x)/w(x)}'. (21) 

In order that V(x) be non-singular, w(x) must be nodeless. Since w(x) is an 
increasing monotonic function (see (19)), the arising of zeros is avoided if 

lim w(x) = or lim w(x) = 0. (22) 

x— >oo x^-oo 

To ensure this requirement it is sufficient that either 

lim uAx) = or lim uAx) = 0. (23) 

x^oo x^-oo 

Such solutions are appropriate for generating non-singular potentials V(x). 
Notice that a similar treatment can be designed for systems defined in a generic 
interval i6(a,i)cl by identifying in (22-23) — oo with a and oo with b. 

4 Complex potentials with real spectrum 



Although in principle irreducible, let us decompose the non-singular 2-SUSY 
transformations of the previous section into two 1-SUSY steps: 



HA 2 = A 2 H 1 , H X A X = AiH, (24) 
where 



d 2 d 
H 1 = -— + V 1 {x), A = — + a l {x), t = l,2. (25) 

The 1-SUSY treatment implies that «i, a 2 obey the Riccati equations: 



-a[ + a 2 ± = V(x) - ei, (26) 

-a , 2 + a$ = V 1 (x)-e 1 , (27) 

where Vi(x) = V(x) + 2a[ and thus V(x) = V(x) + 2(«i + a 2 )' ■ A simple 
comparison of (13) with (26) leads to 

ai(x) = -A (x) = -u'^x) / u^x) , (28) 
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ui(x) being a solution of (18) behaving asymptotically as in (23). Moreover, 
by expanding A = A 2 Ai and comparing the result with (4) we find that: 



a 2 = -a ± + 7] = (3 1 + ( ei - ^/(A - (29) 

This is a particular case of the finite difference Backlund algorithm [3,13,14], 
which algebraically determines a solution to (27) in terms of solutions of (26) 
for two different factorization energies (here e\ and ei). It is interesting as well 
to factorize the involved Hamiltonians: 



H = A[A 1 + e 1 , H 1 = A 1 A^ + e 1 , (30) 
H 1 = A 2 A 2 + e u H = A 2 A- + e u (31) 

A~ = —d/dx + (Xi(x), % — 1,2. Since a±,a 2 and t\ are complex, the A~'s are 
not adjoint to the A£% but A\ = —d/dx + di(x), i — 1,2. 

It is clear now that the complex intermediate potential Vi(x) is non-singular: 



V l (x) = V(x)-2[u' 1 (x)/u l (x)} / . (32) 

To analyse the normalizability of the corresponding eigenfunction associated 
to E n , 



$n( x ) = C nAii) n = C n [lp' n (x) - ufaWnW/u^x)], (33) 

we will employ the operator relationship: 



riA 1 = H-e 1 + A. (34) 

From the validity of (23) and the assumption of \\ip n \\ = 1, it turns out that 
ip n = Aip n /\E n — 6i\ is normalized, and therefore the function r\A\^) n = (E n — 
+ \E n — ei\ip n is normalizable as well. Thus, for A\ty n to be normalizable 
it is neccessary that rf l does not destroy the normalizability of (E n — e^)ty n + 
\E n — €i\ip n . If this is the case (and this will happen for the examples we 
discuss below), we obtain a complex potential V\(x) with real eigenvalues 
E n [21,22,23]. Let us remark that complex Hamiltonians with real spectra have 
been studied recently in the context of PT-symmetry and pseudo-Hermiticity 
[24,25]. 
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5 Illustrative examples 



We shall show that the previous techniques admit very simple applications. 

i) Consider firstly the free particle for which V(x) = 0. The general solution 
Ui(x) of the Schrodinger equation (18) for ei G C is a linear combination of 



e ±{k 1+ ik 2 )x^ ( 35 ) 

where e\ = — (ki+ik 2 ) 2 , k\ > 0, k 2 G IR. In general, such a U\(x) does not tend 
to zero neither when x — > — oo nor when x — > +oo. However, two particular 
solutions with the required behavior are precisely those of (35). We use them 
for obtaining the nodeless w(x): 



w(x) = ±e ±2kix /(2h). (36) 

It turns out that V(x) becomes again the null potential for both w(x), V(x) = 
0. The intermediate 1-SUSY complex potentials generated by using the two 
u\(x) of (35) are as well trivial, V±(x) = 0. Our conclusion is that the null 
potential can be non-trivially transformed in frames of our algorithm only at 
the price of creating singularities (compare with [19]). 

ii) Consider now the well known one-soliton potential (Poschl- Teller) [26] 



V(x) = -2k 2 sech 2 (k x) (37) 

which is obtained from the null potential by a 1-SUSY transformation em- 
ploying cosh(/c x), k > 0. The spectrum of (37) consists of a continuous part 
E > and a bound state at E = —k\ with eigenfunction given by: 



ipo(x) = \J k /2 sech(/c x). (38) 

The general solution ui(x) of (18) for (37) with ei = — (ki+ik 2 ) 2 , ki > 0, k 2 G 
M is a linear combination of the 1-SUSY transformed eigenf unctions of (35) 



e ±ikl+ik2)x [k ttmh(kox) =F (h + ik 2 )}. (39) 

The solutions (39) are precisely the required ones: if the upper signs are taken, 
then ui(x) — * for x — > — oo, while the lower signs ensure Ui(x) — > when 
x — > +oo. An explicit calculation leads to the two nodeless w(x): 
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ke ±2klX cosh[k (x T x )} 

W(x) = ± — r , 40 

V ; 2k x cos\i(k x) ' v ; 

where k$ + k\ + k\ = A; cosh(/c x ), 2/c A;i = k sinh(/c ^o) ; A;, rr £ The two 
2-SUSY partner potentials of (37) read now: 



V(x) = -2klsech 2 [k (x =f x )], (41) 

obtaining just real x -displaced copies of (37). This result has to do with 
the Darboux invariance phenomenon recently discovered for the one-soliton 
well [17,18]. 

On the other hand, the two complex intermediate potentials generated by (39) 
become: 



V 1 (x) = -2klsech 2 [k (x t a>i)], (42) 

where now ki+ik 2 = k cosh(/c a;i), k = k sinh(/c a;i) define 'complex displace- 
ments' Xi G C, kgC. These potentials have a bound state at E = — Ajq whose 
normalized 'ground state' eigenfunction is obtained from (33) by employing 
the ipo(x) of (38) and the u\(x) of (39): 



/ 1 / \ k 
^0 0*0 = ft 



1 /ko + h\ 1 /ko-k^ 
— arctan I — + — arctan ( — 



sech[A;o(x =F xi)](43) 



A; 2 ^ k 2 J k 2 v fc 2 

The complex potentials (42) were obtained for the first time in [22]. 



in) Our final example is the harmonic oscillator: 



V(x) = x 2 , (44) 

which has a purely discrete equidistant spectrum {E n = 2n+ 1, n — 0, 1, . . . }. 
The general solution of (18) for ei G C is now (see [4] and references therein): 



U\{x) = C\t 2 



where iFi(a,c;z) is the Kummer hypergeometric series. In general, such a 
Ui(x) does not satisfy neither lim ;E ^_ 00 Ui(x) = nor lim^oo U\(x) = 0. How- 
ever, there are two particular values for v [y = ±1) leading to solutions with 
the required behavior: 



lfl (i^l,i;^) + 2,S 1 F 1 (i^l,?;^),(45) 
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ui(x) = e 2 



3-ei 3 
~~ 4~ '2 5 



x 



(46) 



Take, e.g., (46) with the upper sign, which in the negative semiaxis x = —\x\ < 
reduces to: 



«iw = T(iy~ e ^^("T" 1 ' ^ ; |:r|2 ) , (47) 

where the Tricomi function \l/(a, c; z) is related with \Fi(a, c; z) through (see, 
e.g., [27]): 



c; *) = F , r(1 f^ i^iK c; z) + ^ ^ z^F^a - c + 1, 2 - c; z) (48) 
1 (a — c + 1) 1 (a) 

Since the leading term in the asymptotic expansion for ty(a,c;z) is z~ a [28], 
one can check that lim^.oo u±(x) = 0. Similarly, if the lower sign is chosen 
we have lim^oo ui(x) = 0. 





20 
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Fig. 1. The real potential V(x) (black curve) generated from V(x) = x 2 (gray curve) 
by means of a 'complex' 2-SUSY transformation with e± = 10 + O.li and the u\{x) 
of (46) with the lower minus sign. 

Once we have identified the solutions (46) with the right asymptotic behavior, 
we evaluated w(x) and then V(x). The resulting expressions in this case are 
too long; instead, we are plotting V(x) for e± — 10 + O.li using (46) with 
the lower minus sign (see figure 1). Contrasting with the results for the pre- 
vious examples, in this case the potentials V(x) are in general different from 
V(x). This means that the transformations involving a pair of complex conju- 
gate factorization energies are effective tools in generating isospectral 2-SUSY 
partner potentials. As a byproduct, we have obtained in a simple way complex 
potentials V\ (x) given by (32) with real energy eigenvalues E n = 2n + l. A plot 
of the 'ground state' probability density |?/>q(x)| 2 , illustrating the existence of 
these bound states for the complex 1-SUSY partner V\ (x) of the oscillator, is 
shown in figure 2. 
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Fig. 2. The 'ground state' probability density |V>q(x)| 2 for the complex 1-SUSY 
partner potential (32) of the oscillator generated by using u\(x) of (46) with the 
lower sign and e± = 10 + O.li 

6 Conclusions 

We have shown that the 2-SUSY transformations involving two complex con- 
jugate factorization energies can produce new non-singular potentials isospec- 
tral to a given initial one. This non-singular character is shared as well by 
the intermediate complex potentials arising when those transformations are 
factorized. 
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